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Abstract :  Recent  data  from  Time  Domain  Pulse 
Ref lectometry  experiments  on  a  three  turn  coil  in  the 
form  of  a  race  track  corroborate  the  theory  of  coupled 
wave  model  for  large  coils.  The  data  demonstrate  ef¬ 
fects  of  turn-to-turn  coupling,  pulse  travel  time,  re¬ 
duced  damping  at  liquid  Nitrogen  temperature,  and  the 
relation  between  short-term  and  long-term  transient  be¬ 
havior.  The  use  of  large  superconductor  coils  as 
energy  storage  elements  inherently  poses  the  problem  of 
dealing  with  fast  electromagnetic  transients  in  the 
vicinity  of  superconductor  boundaries .  Treatment  of 
such  coils  as  the  lumped  inductors  of  circuit  theory  is 
grossly  inaccurate  when  events  to  be  reckoned  are  on 
nanosecond  scale.  Rigorous  treatment  from  the  stand¬ 
point  of  electrodynamic  field  theory  has  resulted  in 
the  Coupled  Wave  Model  which  is  based  on  the  theory  of 
multiple  conductor  transmission  structures.  The  exper¬ 
imental  data  illustrate  the  merit  of  the  model  for 
analysis  of  transients  on  large  coils  and  motivate  fur¬ 
ther  refinement  and  extension. 

Introduction 

One  of  the  important  applications  of  superconduc¬ 
tivity  is  found  in  large  dc  magnets,  which  might  be 
utilized  as  the  primary  source  of  very  strong  steady 
magnetic  fields  or  as  an  energy  storage  element  in 
pulse  power  generators1.  The  large  scale  of  such  mag¬ 
nets  in  all  respects  — physical  size,  current,  and  mag¬ 
netic  energy  — introduce  new  problems  both  in  electro¬ 
dynamic  and  thermal  behavior  which  normally  do  not  pose 
difficulties  in  small  scale  applications,  though  they 
might  be  present  to  a  lesser  degree.  One  such  problem, 
which  the  designer  of  any  large  magnet  must  face,  is 
electromagnetic  transients  and  the  likelihood  of  high 
voltage  breakdown  together  with  the  damage  that  results 
from  them. 

The  occurrence  of  high  voltage  breakdown  between 
layers  of  windings  of  large  coils  is  not  new  since  it 
has  been  one  of  the  lingering  problems  in  the  design  of 
large  power  transformers.  That  breakdown  should  occur 
in  alternating  current  transformers  is  not  so  surpris¬ 
ing  from  a  theoretical  standpoint.  However,  at  first 
glance,  it  would  seem  very  unusual  that  it  should  also 
occur  in  direct  current  superconducting  magnets.  Yet, 
evidence  of  high  voltage  breakdown  has  been  observed  in 
such  magnets  following  emergency  shut  down  when  the 
stored  magnetic  energy  is  dumped  into  shunt  resistors 
in  accordance  with  the  conventional  way  of  limiting 
voltages  across  inductors.  In  the  conventional  view¬ 
point,  a  coil  of  wire  is  treated  as  the  lumped  inductor 
of  circuit  theory.  According  to  this  theory,  during 
the  transient  current  decay  period,  the  voltage  across 
the  terminals  and  anywhere  within  the  coil  should  never 
exceed  the  voltage  across  the  limiting  shunt  resistor. 
Thus,  it  is  clear  that  the  conventional  circuit  repre¬ 
sentation  of  a  coil  in  not  adequate  to  account  for  high 
voltages.  A  better  understanding  of  electromagnetic 
transients  is  needed. 

The  inadequacy  of  the  conventional  approach  lies 
in  the  fact  that  any  circuit  representation  of  a  phys¬ 
ical  system,  by  nature  having  no  provision  for  spatial 
distribution,  treats  all  events  simultaneously.  In 
reality,  however,  electromagnetic  events  do  not  occur 


simultaneously  at  any  pair  of  points  separated  by  a 
distance.  The  bearing  of  this  fact  on  the  transient 
behavior  of  large  coils  is  appreciated  when  one  consi¬ 
ders  that  an  electromagnetic  disturbance  travels  one 
foot  in  one  nanosecond.  Consequently,  circuit  theory 
is  grossly  inappropriate  to  attack  transients  on  large 
scale.  It  soon  became  apparent  that  the  key  to  under¬ 
standing  such  transients  lies  in  field  theoretic  model¬ 
ling.  In  an  earlier  paper,  Gabriel  and  Burkhart2  re¬ 
ported  observation  of  multiplicity  of  resonances  in 
frequency  as  well  as  evidence  of  multiple  wave  reflec¬ 
tions  in  time.  A  very  crude  analysis  based  on  treating 
the  coil  as  a  single  wire  transmission  line  yielded 
surprisingly  good  qualitative  agreement.  Based  on  this 
success,  attention  was  directed  to  laying  the  founda¬ 
tions  for  the  more  realistic  Coupled  Wave  Model  the  de¬ 
tails  of  which  are  given  in  an  earlier  document3.  In 
this  paper,  a  brief  account  is  given  of  further  devel¬ 
opment  of  the  model,  and  recent  data  are  presented  as  a 
preliminary  step  toward  corroboration  and  further  evo¬ 
lution. 

Theoretical  Foundation 

Outline  of  Model 

The  objective  has  been  to  evolve  a  model  which  is 
both  rigorously  founded  on  field  theory  and  broadly 
compatible  with  experimental  evidence.  The  most  sug¬ 
gestive  evidence  is  the  observed  travel  time  of  a  dis¬ 
turbance  around  a  single  turn,  because  it  indicates 
that  the  transient  field  behaves  very  nearly  as  a  Tran- 
verse  ElectroMagnetic  (TEM)  wave.  Motivated  by  this 
fact,  the  model  is  founded  on  the  theory  of  multiple 
conductor  transmission  lines  which  support  TEM  waves, 
even  though  the  geometry  of  any  coil  is  not  exactly 
compatible  with  the  conditions  for  TEM  waves  In  the 
strict  sense.  The  gist  of  the  model  is  conveyed  by  a 
wave  coupling  matrix,  conveniently  called  an  inductance 
or  capacitance  matrix,  and  the  crucial  interconnection 
matrix  containing  the  naturally  periodic  boundary  con¬ 
ditions  demanded  by  a  coil  structure.  With  C  measuring 
the  coordinate  distance  along  any  single  turn  of  length 
2L,  the  electric  and  magnetic  fields  on  the  surface  of 
the  wire  of  the  i*1*1  turn  are  characterized  by  amplitude 
functions  Vf(C,t)  and  ii(C,t),  conveniently  called  vol¬ 
tage  and  current.  A  coil  of  M  turns  is  then  described 
by  the  M-dimensional  column  vectors  v(C»t)  and  i(C,t) 
which  satisfy  the  standard  telegraphist  equations 
9  9 

—  v(C,t)  =  -lay"1  — i(C,t)  (1) 

3C  3t 

9  9 

-i(!,t)=-n  —  v(c,t)  (2) 

_  3c  3t 

where  y  is  the  symmetric  coupling  matrix.  These  equa¬ 
tions  result  from  appropriate  normalization  of  the 
field  equations  for  the  vector  potential  A(C,t)  and  the 
scalar  Coulomb  potential  <KC,t),  which  are  actually 
components  of  the  Lorentz  invariant  potential  in  4- 
space . 

After  Laplace  transformation,  and  denoting  the 
transforms  as  V(C,s)  and  I(C,s),  one  finds  the  general 
solution  to  be,  as  expected,  a  set  of  waves  travelling 
in  unison,  viz. 
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where  A  and  15  are  column  vectors  to  be  determined  by 
boundary  conditions,  and  Y  =  /e/p  y  is  a  characteristic 
wave  admittance  which  is  proportional  to  the  geometric 
coupling  matrix  through  the  permeability  v  and  permiti- 
vitty  e  of  the  surrounding  medium. 

So  far  the  development  is  quite  standard  and  well 
known.  However,  what  distinguishes  the  coupled  wave 
model  from  the  ordinary  multiple  conductor  transmission 
line  is  the  boundary  conditions  and  their  representa¬ 
tion  by  an  interconnection  matrix  T.  In  a  uniformly 
wound  coll  of  M  turns,  regardless  of  shape,  the  physi¬ 
cal  input  terminals  of  the  coil  as  a  unit  are  taken  to 
be  at  ?=-L  on  the  first  turn  and  at  ?=L  on  the  M11*1 
turn,  Vt  and  It  being  the  terminal  voltage  and  current 
demanded  by  external  connections .  The  periodic  struc¬ 
ture  of  the  coil,  on  the  other  hand,  demands  equality 
of  the  voltage  and  current  at  ?=-L  on  the  m6*1  turn  to 
those  at  C=L  on  the  (m-l)*^  turn.  In  brief,  these  con¬ 
ditions  are  effectively  expressed  by  the  vector  equa¬ 
tions 

T  I(-L,s)  -  1(1,, s)  =  0  (5) 

T  V(-L,s)  -  V(L,s)  =  VtaM  (6) 


the  coil  were  viewed  as  a  lumped  inductance  L0,  its  ad¬ 
mittance  would  be  l/sL0.  Even  a  representation  by  a 
finite  number  of  hypothetical  capacitance-inductance 
meshes,  as  commonly  employed  in  analysis  of  transfor¬ 
mers,  would  still  yield  an  admittance  which  is  the 
ratio  of  rational  polynomials  in  s .  By  constrast,  the 
admittance  predicted  by  the  coupled  wave  model  is  a 
function  not  of  s  but  of  the  delay  parameter  e-sT, 
where  x  =  2L/c  is  the  travel  time  around  a  single  turn. 
It  can  be  shown  that  in  general  the  admittance  may  be 
expanded  in  a  series 

00 

Yt  =  a0  +  I  an  e-"sT  (8) 

n=l 

where  an  are  constant  coefficients  determined  by  the 
coupling  matrix  elements  Y-tj.  Significantly,  the  expo¬ 
nential  terms  indicate  that  voltages  and  currents,  as 
observables,  evolve  not  as  continuous  functions  in  time 
but  in  discrete  steps  having  durations  equal  to  the 
characteristic  travel  time  around  a  single  turn.  Even 
in  the  case  of  three  turns,  however,  calculation  of  the 
coefficients  an  for  a  given  matrix  Y  is  not  a  simple 
task.  Nevertheless,  the  first  four  coefficients  have 
been  derived  for  the  purpose  of  comparing  to  the  exper¬ 
imental  observations  in  broad  terms.  In  this  case,  the 
coil  is  appropriately  characterized  by  the  characteris¬ 
tic  wave  admittance  matrix  which  takes  the  general 
form 


where 


Y  = 


Y1+Y3 

-Y2 

-y3 


-y2  — y3 

Yi+Y2  .-y2 
— Y2  Y1+Y3 


(9) 


*M  = 


0 

1 


and  T  is  the  interconnection  matrix.  In  a  uniformly 
wound  coil,  with  no  folding  layers,  it  takes  the  parti¬ 
cular  form 
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while  the  interconnection  matrix  reduces  to 


0 

0 

1 


1  0 
0  1 
0  0 


(10) 


The  3-turn  coil  is  useful  for  the  purpose  of  prelimin¬ 
ary  experimental  verification  because  it  does  not  suf¬ 
fer  from  degeneracy  as  the  2-turn  coil  while  at  the 
same  time  it  holds  the  calculations  to  reasonable  pro¬ 
portions  . 


Theory  of  Experimental  Method 


More  generally,  every  uniform  coil,  including  those 
with  folding  layers,  has  an  interconnection  matrix 
whose  form  is  such  that  every  column  and  every  row  has 
but  one  nonzero  element.  Such  matrices  are  orthogonal 
having  eigenvalues  that  lie  on  the  unit  circle  in  the 
complex  plane  which  are  intimately  related  to  the  re¬ 
sonances  normally  observed  in  any  coil. 

The  two-turn  coil  is  the  simplest  symmetric  con¬ 
figuration  which  permits  straight  forward  exact  solu¬ 
tion  of  the  boundary  value  problem  consisting  of  Eqs. 
(l)-(6).  Such  solution  is  analyzed  in  detail  in  the 
earlier  work3.  The  general  case  of  M-turn  coil,  how¬ 
ever,  poses  formidable  difficulties  owing  to  the  neces¬ 
sity  for  tedious  algebraic  manipulation  of  MxM  ma¬ 
trices.  Nevertheless,  some  broad  features  can  be  de¬ 
duced  without  explicit  calculation. 

A  significant  parameter,  especially  for  our  pur¬ 
pose  here  in  light  of  the  experimental  data,  is  the  in¬ 
put  terminal  admittance  Yt  since  it  plays  a  dominant 
role  in  the  way  the  coil  is  interfaced  with  external 
systems.  As  usual,  it  is  defined  as 

Yt  -  it  .  .  (7) 

Vt  Vt(s) 


Observation  of  the  effects  of  the  terminal  admit¬ 
tance  function  offers  an  excellent  route  to  verifica¬ 
tion  and  appraisal.  Thus,  the  core  of  the  experimental 
method  consists  of  time  domain  ref lectrometry  (TDR) 
both  in  pulse  and  step  form.  It  is  well  known  that 
when  a  cable  is  terminated  with  an  admittance  Yp(s)  as 
the  load,  the  reflection  coefficient  in  terms  of  nor¬ 
malized  admittance  is  given  by 

P(s)  =  —  ,  yt  -  Yt/Y0  (11) 

l+yt 

where  Y0  is  the  characteristic  admittance  of  the  cable, 
in  this  case  1/50  mho.  Moreover,  if  Vg(s)  is  the  equi¬ 
valent  voltage  (as  amplitude  of  the  electric  field  at  a 
reference  terminal  plane)  of  a  matched  generator,  then 
the  reflected  voltage  observed  at  the  generator  ter¬ 
minal  is  given  by 

Vr(s)  =  p(s)Vg(s).  (13) 

Obviously,  in  time  domain  vr(t)  is  a  convolution  of 
Vg(t)  with  the  inverse  of  P(s).  Alternately,  the  in¬ 
version  can  be  effected  by  using  the  series  form  of  the 
admittance  function  (8)  to  expand  the  reflection  coef¬ 
ficient  p(s)  in  (11)  into  another  series 

00 

p(s)  -  p0[l  +  I  Pne-nST]  • 
n=l 


The  remarkable  departure  of  the  coupled  wave  model  from 
circuit  theory  appears  in  this  admittance  function.  If 
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(14) 


If  Vg(t)  is  arranged  to  be  non-zero  only  over  an  inter¬ 
val  shorter  than  t,  then  a  term-by-term  inversion  of 
(13)  using  (14)  yields  the  sequence 

00 

vr(t)  =  P0Vg(t)  +  0o  I  PnVg(t-nt).  (15) 

n-1 

For  the  case  Vg(t)=u(t),  the  unit  step  function,  the 
inversion  becomes  a  sequence  of  weighted  but  succes¬ 
sively  delayed  steps 

00 

vr(t)  “  P0u(t)  +  P0  I  Pn  u(t-nT).  (16) 
n*l 

Equations  (15)  and  (16)  are  the  analytical  expres¬ 
sions  of  the  primary  observables  presented  below,  and 
they  emphasize  the  significant  departure  from  the  cir¬ 
cuit  theoretic  inductance  model.  In  particular,  if  the 
circuit  admittance  Yt=l/sL0  were  employed  in  (11)  and 
(13),  the  reflected  voltage  for  unit  step  input  would 
be  a  single  decaying  exponential  which  is  far  different 
from  the  sequence  (16)  on  short  time  scale.  It  is  re¬ 
markable,  however,  to  note  that  because  pn  are  monoton- 
ically  decreasing,  on  a  coarse  time  scale  larger  chan 
T,  the  sum  In  (16)  approaches  an  exponential  decay 
whose  time  constant  is  determined  by  the  system  resis¬ 
tance  which  in  this  case  Is  principally  the  generator 
resistance  rather  than  the  much  smaller  coil  resis¬ 
tance.  This  asymptotic  behavior  has  been  verified  by 
computer  simulation  and  is  born  out  in  the  data  below. 
For  the  sake  of  completeness,  the  first  four  coeffi¬ 
cients  which  are  calculated  from  the  matrices  (9)  and 
(10)  for  the  3-turn  coil  are  listed  here  in  terms  of 
normalized  elements  yi  =  Yi/Y0: 


Po  =  (2-y0)/(2+y0) 

Pi  * 

8y2/p 

2 

p2  =  4y2  y3/p 

P3  “ 

-8y0/p 

2 

Yo  =  Yl  +  Y3 

P  = 

4  -  y 

o 


Note  that  the  second  and  third  terms,  p^  and  P2,  are 
positive  and  depend  on  the  mutual  turn-to-turn  coupling 
while  the  fourth  term  is  negative  and  principally  de¬ 
pendent  on  self  coupling  corresponding  to  the  diagonal 
terms  of  the  Y  matrix. 

TDR  Experiments 

Time  domain  reflectrometry  experiments  were  con¬ 
ducted  on  a  3-tum  coil  consisting  of  ordinary  No.  22 
enameled  copper  wire.  The  coil  is  in  the  shape  of  a 
race  track  (or  paper  clip)  lying  in  a  plane  and  measur¬ 
ing  100  cm  in  length  and  15  cm  in  width,  the  ends  being 
semicircles  with  15  cm  diameter.  The  overall  peripher¬ 
al  length  of  each  turn  amounts  to  220  cm  which  is  the 


equivalent  of  approximately  7.3  nanoseconds  travel  time 
per  turn.  Choice  of  these  dimensions  is  dictated  by 
considerations  of  pulse  resolution,  given  that  the 
pulse  width  of  the  generator  used  is  approximately  5 
nanosecond. 

A  schematic  diagram  for  the  TDR  measurements  is 
shown  in  Figure  1.  Briefly,  a  voltage  pulse  emerging 
from  the  generator  is  split  in  two,  one  half  going  to 
the  coil  along  the  main  line,  and  the  other  half  going 
to  the  oscilloscope  along  the  auxiliarly  line.  The 
fixed  delay  suffered  along  these  lines  is  common  to  all 
pulses  and  thus  it  is  inconsequential.  The  main  re¬ 
sults  are  summarized  next: 

1.  In  Figures  2a  and  2b,  the  leftmost  faint  trace 
Is  the  incident  pulse.  The  next  pulse  from  the  left  is 
the  first  reflected  pulse  due  to  the  discontinuity  at 
the  coil  terminals,  and  it  is  proportional  to  the  first 
term  p0  in  Eq.  (15).  If  there  were  no  tum-to-turn 
coupling,  the  incident  pulse  entering  terminal  1  as 
positive  and  terminal  2  as  negative  travels  the  entire 
course  of  three  turns  amounting  to  a  delay  of  approxi¬ 
mately  22  nanosecond,  arriving  back  at  the  terminals  In 
inverted  polarity,  and  thence  back  to  the  generator  and 
oscilloscope.  This  appears  as  the  inverted  pulse  in 
the  trace,  and  it  corresponds  to  the  fourth  term  P3  in 
Eq.  (15).  Because  of  coupling,  however.  Immediately 
upon  arrival  at  the  outer  turns,  the  incident  pulse 
induces  secondary  ones  having  reversed  polarity  in  the 
inner  second  and  third  turns .  These  secondary  pulses 
travel  along  one  turn  and  two  turns,  respectively,  be¬ 
fore  arriving  at  the  terminals  on  their  way  towards  the 
generator  in  the  same  polarity  as  the  Incident  pulse. 
The  result  is  the  two  positive  echo  pulses  appearing  in 
the  traces  at  approximately  7.5  and  15  nanosecond  after 
the  first  reflection.  These  correspond  to  the  terms  Pi 
and  P2  in  Eq.  (15)  which  are  dependent  on  the  off- 
diagonal  elements  of  the  coupling  matrix.  It  is  inter¬ 
esting  to  note  that  coupling  beyond  nearest  neighbors 
is  significant,  because  if  Y3  and  therefore  y3  were 
zero  the  echo  at  15  nanosecond  would  not  be  present. 

2.  The  step  response  on  short  time  scale  appears 
In  Figure  3a  and  3b.  The  events  are  similar  to  the 
pulse  case,  save  for  the  fact  that  the  sequence  is 
cumulatively  additive.  The  magnitude  of  each  of  the 
stpes  is  proportional  to  the  respective  terms  pn  in  Eq. 
(16). 

3.  Figures  4a  and  4b  illustrate  in  a  spectacular 
way  the  asymptotic  behavior  of  the  coil  over  long  time 
scale,  which  approaches  the  exponential  decay  expected 
from  circuit  theory.  In  the  wave  viewpoint,  however, 
the  rate  of  decay  is  essentially  a  manifestation  of  the 
rate  of  decrease  of  the  coefficients  Pn  as  n  increases 
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Figure  1:  Schematic  Diagram  of  TDR  Arrangement  and  Coil 


612 


Figure  3:  Step  reflection,  a)  Room  temperature;  b)  Liquid  nitrogen  temperature 


(a)  (b) 

Figure  4:  Asymptotic  long  time  reflection,  a)  Room  temperature;  b)  Liquid  nitrogen  temperature 
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indefinitely.  This  rate  is  determined  by  the  ratio  of 
the  admittance  elements  to  the  generator  conductance 
(which  is  matched  to  the  cable  admittance  Y0).  If  the 
coil  and  generator  are  treated  from  a  circuit  stand¬ 
point,  taking  L0  in  series  with  the  generator  resis¬ 
tance  l/YOJ  then  clearly  the  decay  time  constant  would 
indeed  be  L0Y0. 

4.  Operation  at  liquid  Nitrogen  temperature  re¬ 
duces  the  finite  resistance  of  the  coil  wire.  The  re¬ 
sulting  reduction  in  damping  is  mainly  evident  in  the 
short  time  behavior  of  the  coil  as  shown  in  Figure  2b 
for  pulse  input,  and  in  Figure  3b  for  step  input.  The 
effect  is  most  pronounced  in  the  inverted  echo  which, 
as  indicated  above,  travels  the  entire  length  of  3- 
turns.  This  echo  is  principally  determined  by  the  ele¬ 
ment  yQ  of  the  admittance  matrix.  Theoretical  consid¬ 
erations  based  on  a  minimum  energy  principle1*  indicate 
that  the  effect  of  finite  wire  resistance  should  appear 
as  a  first  order  additive  correction  on  only  the  diag¬ 
onal  elements  of  the  matrix  which  is  consistent  with 
the  observed  result.  On  the  other  hand,  Figures  4a  and 
4b  show  no  discernible  difference  in  the  long  term  de¬ 
cay  rates  at  the  two  different  temperatures.  As  indi¬ 
cated  in  3  above,  the  generator  resistance  of  50  ohm, 
being  much  larger  than  the  coil  resistance  both  at  room 
and  at  liquid  Nitrogen  temperatures,  is  the  dominant 
resistance  in  determining  the  long  term  decay  rate. 
Thus,  this  result  is  also  consistent. 

Conclusion 

In  summary,  the  travel  time  of  a  disturbance 
around  a  single  turn  and  the  turn-to-turn  coupling  are 
the  dominant  aspects  by  which  the  Coupled  Wave  Model 
departs  from  circuit  theory  in  describing  transient  re¬ 
sponse  of  coils  over  short  time  scale.  The  excellent 
agreement  of  the  model  with  the  experimental  observa¬ 
tions  enumerated  above  motivates  confidence  in  the 


basic  suppositions  undelying  its  theory.  As  to  the 
problem  of  high  voltage  breakdown,  it  is  most  likely  to 
occur  on  the  short  time  scale  comparable  to  travel  time 
per  turn,  as  indicated  by  computer  calculation  on  a 
two-turn  coil3.  Success  of  the  model  thus  far  makes  it 
well  worth  the  effort  to  tackle  the  more  complicated 
and  difficult  configurations  with  a  view  toward  devel¬ 
oping  refinements  that  would  permit  prediction  and  sim¬ 
ulation. 
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